A Principal Components Analysis and Functional Annotation of Differentially Expressed Genes in Brain Regions of Gray Rats Selected for Tame or Aggressive Behavior.

The process of domestication, despite its short duration as it compared with the time scale of the natural evolutionary process, has caused rapid and substantial changes in the phenotype of domestic animal species. Nonetheless, the genetic mechanisms underlying these changes remain poorly understood. The present study deals with an analysis of the transcriptomes from four brain regions of gray rats (Rattus norvegicus), serving as an experimental model object of domestication. We compared gene expression profiles in the hypothalamus, hippocampus, periaqueductal gray matter, and the midbrain tegmental region between tame domesticated and aggressive gray rats and revealed subdivisions of differentially expressed genes by principal components analysis that explain the main part of differentially gene expression variance. Functional analysis (in the DAVID (Database for Annotation, Visualization and Integrated Discovery) Bioinformatics Resources database) of the differentially expressed genes allowed us to identify and describe the key biological processes that can participate in the formation of the different behavioral patterns seen in the two groups of gray rats. Using the STRING- DB (search tool for recurring instances of neighboring genes) web service, we built a gene association network. The genes engaged in broad network interactions have been identified. Our study offers data on the genes whose expression levels change in response to artificial selection for behavior during animal domestication.

Age-Dependent Changes in the Relationships between Traits Associated with the Pathogenesis of Stress-Sensitive Hypertension in ISIAH Rats.

Hypertension is one of the most significant risk factors for many cardiovascular diseases. At different stages of hypertension development, various pathophysiological processes can play a key role in the manifestation of the hypertensive phenotype and of comorbid conditions. Accordingly, it is thought that when diagnosing and choosing a strategy for treating hypertension, it is necessary to take into account age, the stage of disorder development, comorbidities, and effects of emotional-psychosocial factors. Nonetheless, such an approach to choosing a treatment strategy is hampered by incomplete knowledge about details of age-related associations between the numerous features that may contribute to the manifestation of the hypertensive phenotype. Here, we used two groups of male F2(ISIAHxWAG) hybrids of different ages, obtained by crossing hypertensive ISIAH rats (simulating stress-sensitive arterial hypertension) and normotensive WAG rats. By principal component analysis, the relationships among 21 morphological, physiological, and behavioral traits were examined. It was shown that the development of stress-sensitive hypertension in ISIAH rats is accompanied not only by an age-dependent (FDR < 5%) persistent increase in basal blood pressure but also by a decrease in the response to stress and by an increase in anxiety. The plasma corticosterone concentration at rest and its increase during short-term restraint stress in a group of young rats did not have a straightforward relationship with the other analyzed traits. Nonetheless, in older animals, such associations were found. Thus, the study revealed age-dependent relationships between the key features that determine hypertension manifestation in ISIAH rats. Our results may be useful for designing therapeutic strategies against stress-sensitive hypertension, taking into account the patients' age.

Candidate SNP Markers Significantly Altering the Affinity of TATA-Binding Protein for the Promoters of Human Hub Genes for Atherogenesis, Atherosclerosis and Atheroprotection.

Atherosclerosis is a systemic disease in which focal lesions in arteries promote the build-up of lipoproteins and cholesterol they are transporting. The development of atheroma (atherogenesis) narrows blood vessels, reduces the blood supply and leads to cardiovascular diseases. According to the World Health Organization (WHO), cardiovascular diseases are the leading cause of death, which has been especially boosted since the COVID-19 pandemic. There is a variety of contributors to atherosclerosis, including lifestyle factors and genetic predisposition. Antioxidant diets and recreational exercises act as atheroprotectors and can retard atherogenesis. The search for molecular markers of atherogenesis and atheroprotection for predictive, preventive and personalized medicine appears to be the most promising direction for the study of atherosclerosis. In this work, we have analyzed 1068 human genes associated with atherogenesis, atherosclerosis and atheroprotection. The hub genes regulating these processes have been found to be the most ancient. In silico analysis of all 5112 SNPs in their promoters has revealed 330 candidate SNP markers, which statistically significantly change the affinity of the TATA-binding protein (TBP) for these promoters. These molecular markers have made us confident that natural selection acts against underexpression of the hub genes for atherogenesis, atherosclerosis and atheroprotection. At the same time, upregulation of the one for atheroprotection promotes human health.

Differentially Expressed Genes and Molecular Susceptibility to Human Age-Related Diseases.

Mainstream transcriptome profiling of susceptibility versus resistance to age-related diseases (ARDs) is focused on differentially expressed genes (DEGs) specific to gender, age, and pathogeneses. This approach fits in well with predictive, preventive, personalized, participatory medicine and helps understand how, why, when, and what ARDs one can develop depending on their genetic background. Within this mainstream paradigm, we wanted to find out whether the known ARD-linked DEGs available in PubMed can reveal a molecular marker that will serve the purpose in anyone's any tissue at any time. We sequenced the periaqueductal gray (PAG) transcriptome of tame versus aggressive rats, identified rat-behavior-related DEGs, and compared them with their known homologous animal ARD-linked DEGs. This analysis yielded statistically significant correlations between behavior-related and ARD-susceptibility-related fold changes (log2 values) in the expression of these DEG homologs. We found principal components, PC1 and PC2, corresponding to the half-sum and the half-difference of these log2 values, respectively. With the DEGs linked to ARD susceptibility and ARD resistance in humans used as controls, we verified these principal components. This yielded only one statistically significant common molecular marker for ARDs: an excess of Fcγ receptor IIb suppressing immune cell hyperactivation.

The complex role of transcription factor GAGA in germline death during Drosophila spermatogenesis: transcriptomic and bioinformatic analyses.

The GAGA protein (also known as GAF) is a transcription factor encoded by the Trl gene in D. melanogaster. GAGA is involved in the regulation of transcription of many genes at all stages of fly development and life. Recently, we investigated the participation of GAGA in spermatogenesis and discovered that Trl mutants experience massive degradation of germline cells in the testes. Trl underexpression induces autophagic death of spermatocytes, thereby leading to reduced testis size. Here, we aimed to determine the role of the transcription factor GAGA in the regulation of ectopic germline cell death. We investigated how Trl underexpression affects gene expression in the testes. We identified 15,993 genes in three biological replicates of our RNA-seq analysis and compared transcript levels between hypomorphic Trl R85/Trl 362 and Oregon testes. A total of 2,437 differentially expressed genes were found, including 1,686 upregulated and 751 downregulated genes. At the transcriptional level, we detected the development of cellular stress in the Trl-mutant testes: downregulation of the genes normally expressed in the testes (indicating slowed or abrogated spermatocyte differentiation) and increased expression of metabolic and proteolysis-related genes, including stress response long noncoding RNAs. Nonetheless, in the Flybase Gene Ontology lists of genes related to cell death, autophagy, or stress, there was no enrichment with GAGA-binding sites. Furthermore, we did not identify any specific GAGA-dependent cell death pathway that could regulate spermatocyte death. Thus, our data suggest that GAGA deficiency in male germline cells leads to an imbalance of metabolic processes, impaired mitochondrial function, and cell death due to cellular stress.

Transcription Factors as Important Regulators of Changes in Behavior through Domestication of Gray Rats: Quantitative Data from RNA Sequencing.

Studies on hereditary fixation of the tame-behavior phenotype during animal domestication remain relevant and important because they are of both basic research and applied significance. In model animals, gray rats Rattus norvegicus bred for either an enhancement or reduction in defensive response to humans, for the first time, we used high-throughput RNA sequencing to investigate differential expression of genes in tissue samples from the tegmental region of the midbrain in 2-month-old rats showing either tame or aggressive behavior. A total of 42 differentially expressed genes (DEGs; adjusted p-value < 0.01 and fold-change > 2) were identified, with 20 upregulated and 22 downregulated genes in the tissue samples from tame rats compared with aggressive rats. Among them, three genes encoding transcription factors (TFs) were detected: Ascl3 was upregulated, whereas Fos and Fosb were downregulated in tissue samples from the brains of tame rats brain. Other DEGs were annotated as associated with extracellular matrix components, transporter proteins, the neurotransmitter system, signaling molecules, and immune system proteins. We believe that these DEGs encode proteins that constitute a multifactorial system determining the behavior for which the rats have been artificially selected. We demonstrated that several structural subtypes of E-box motifs­known as binding sites for many developmental TFs of the bHLH class, including the ASCL subfamily of TFs­are enriched in the set of promoters of the DEGs downregulated in the tissue samples of tame rats'. Because ASCL3 may act as a repressor on target genes of other developmental TFs of the bHLH class, we hypothesize that the expression of TF gene Ascl3 in tame rats indicates longer neurogenesis (as compared to aggressive rats), which is a sign of neoteny and domestication. Thus, our domestication model shows a new function of TF ASCL3: it may play the most important role in behavioral changes in animals.

Web-MCOT Server for Motif Co-Occurrence Search in ChIP-Seq Data.

(1) Background: The widespread application of ChIP-seq technology requires annotation of cis-regulatory modules through the search of co-occurred motifs. (2) Methods: We present the web server Motifs Co-Occurrence Tool (Web-MCOT) that for a single ChIP-seq dataset detects the composite elements (CEs) or overrepresented homo- and heterotypic pairs of motifs with spacers and overlaps, with any mutual orientations, uncovering various similarities to recognition models within pairs of motifs. The first (Anchor) motif in CEs respects the target transcription factor of the ChIP-seq experiment, while the second one (Partner) can be defined either by a user or a public library of Partner motifs being processed. (3) Results: Web-MCOT computes the significances of CEs without reference to motif conservation and those with more conserved Partner and Anchor motifs. Graphic results show histograms of CE abundance depending on orientations of motifs, overlap and spacer lengths; logos of the most common CE structural types with an overlap of motifs, and heatmaps depicting the abundance of CEs with one motif possessing higher conservation than another. (4) Conclusions: Novel capacities of Web-MCOT allow retrieving from a single ChIP-seq dataset with maximal information on the co-occurrence of motifs and potentiates planning of next ChIP-seq experiments.

Stress Reactivity, Susceptibility to Hypertension, and Differential Expression of Genes in Hypertensive Compared to Normotensive Patients.

Although half of hypertensive patients have hypertensive parents, known hypertension-related human loci identified by genome-wide analysis explain only 3% of hypertension heredity. Therefore, mainstream transcriptome profiling of hypertensive subjects addresses differentially expressed genes (DEGs) specific to gender, age, and comorbidities in accordance with predictive preventive personalized participatory medicine treating patients according to their symptoms, individual lifestyle, and genetic background. Within this mainstream paradigm, here, we determined whether, among the known hypertension-related DEGs that we could find, there is any genome-wide hypertension theranostic molecular marker applicable to everyone, everywhere, anytime. Therefore, we sequenced the hippocampal transcriptome of tame and aggressive rats, corresponding to low and high stress reactivity, an increase of which raises hypertensive risk; we identified stress-reactivity-related rat DEGs and compared them with their known homologous hypertension-related animal DEGs. This yielded significant correlations between stress reactivity-related and hypertension-related fold changes (log2 values) of these DEG homologs. We found principal components, PC1 and PC2, corresponding to a half-difference and half-sum of these log2 values. Using the DEGs of hypertensive versus normotensive patients (as the control), we verified the correlations and principal components. This analysis highlighted downregulation of ß-protocadherins and hemoglobin as whole-genome hypertension theranostic molecular markers associated with a wide vascular inner diameter and low blood viscosity, respectively.

E. coli aggregation and impaired cell division after terahertz irradiation.

In this study we demonstrated that exposure of Escherichia coli (E. coli) to terahertz (THz) radiation resulted in a change in the activities of the tdcABCDEFGR and matA-F genes (signs of cell aggregation), gene yjjQ (signs of suppression of cell motility), dicABCF, FtsZ, and minCDE genes (signs of suppression of cell division), sfmACDHF genes (signs of adhesin synthesis), yjbEFGH and gfcA genes (signs of cell envelope stabilization). Moreover, THz radiation induced E. coli csg operon genes of amyloid biosynthesis. Electron microscopy revealed that the irradiated bacteria underwent increased aggregation; 20% of them formed bundle-like structures consisting of two to four pili clumped together. This could be the result of changes in the adhesive properties of the pili. We also found aberrations in cell wall structure in the middle part of the bacterial cell; these aberrations impaired the cell at the initial stages of division and resulted in accumulation of long rod-like cells. Overall, THz radiation was shown to have adverse effects on bacterial populations resulting in cells with abnormal morphology.

Domestication Explains Two-Thirds of Differential-Gene-Expression Variance between Domestic and Wild Animals; The Remaining One-Third Reflects Intraspecific and Interspecific Variation.

Belyaev's concept of destabilizing selection during domestication was a major achievement in the XX century. Its practical value has been realized in commercial colors of the domesticated fox that never occur in the wild and has been confirmed in a wide variety of pet breeds. Many human disease models involving animals allow to test drugs before human testing. Perhaps this is why investigators doing transcriptomic profiling of domestic versus wild animals have searched for breed-specific patterns. Here we sequenced hypothalamic transcriptomes of tame and aggressive rats, identified their differentially expressed genes (DEGs), and, for the first time, applied principal component analysis to compare them with all the known DEGs of domestic versus wild animals that we could find. Two principal components, PC1 and PC2, respectively explained 67% and 33% of differential-gene-expression variance (hereinafter: log2 value) between domestic and wild animals. PC1 corresponded to multiple orthologous DEGs supported by homologs; these DEGs kept the log2 value sign from species to species and from tissue to tissue (i.e., a common domestication pattern). PC2 represented stand-alone homologous DEG pairs reversing the log2 value sign from one species to another and from tissue to tissue (i.e., representing intraspecific and interspecific variation).

Disruptive Selection of Human Immunostimulatory and Immunosuppressive Genes Both Provokes and Prevents Rheumatoid Arthritis, Respectively, as a Self-Domestication Syndrome.

Using our previously published Web service SNP_TATA_Comparator, we conducted a genome-wide study of single-nucleotide polymorphisms (SNPs) within core promoters of 68 human rheumatoid arthritis (RA)-related genes. Using 603 SNPs within 25 genes clinically associated with RA-comorbid disorders, we predicted 84 and 70 candidate SNP markers for overexpression and underexpression of these genes, respectively, among which 58 and 96 candidate SNP markers, respectively, can relieve and worsen RA as if there is a neutral drift toward susceptibility to RA. Similarly, we predicted natural selection toward susceptibility to RA for 8 immunostimulatory genes (e.g., IL9R) and 10 genes most often associated with RA (e.g., NPY). On the contrary, using 25 immunosuppressive genes, we predicted 70 and 109 candidate SNP markers aggravating and relieving RA, respectively (e.g., IL1R2 and TGFB2), suggesting that natural selection can simultaneously additionally yield resistance to RA. We concluded that disruptive natural selection of human immunostimulatory and immunosuppressive genes is concurrently elevating and reducing the risk of RA, respectively. So, we hypothesize that RA in human could be a self-domestication syndrome referring to evolution patterns in domestic animals. We tested this hypothesis by means of public RNA-Seq data on 1740 differentially expressed genes (DEGs) of pets vs. wild animals (e.g., dogs vs. wolves). The number of DEGs in the domestic animals corresponding to worsened RA condition in humans was significantly larger than that in the related wild animals (10 vs. 3). Moreover, much less DEGs in the domestic animals were accordant to relieved RA condition in humans than those in the wild animals (1 vs. 8 genes). This indicates that the anthropogenic environment, in contrast to a natural one, affects gene expression across the whole genome (e.g., immunostimulatory and immunosuppressive genes) in a manner that likely contributes to RA. The difference in gene numbers is statistically significant as confirmed by binomial distribution (p < 0.01), Pearson's χ2 (p < 0.01), and Fisher's exact test (p < 0.05). This allows us to propose RA as a candidate symptom within a self-domestication syndrome. Such syndrome might be considered as a human's payment with health for the benefits received during evolution.

A Bioinformatics Model of Human Diseases on the Basis of Differentially Expressed Genes (of Domestic Versus Wild Animals) That Are Orthologs of Human Genes Associated with Reproductive-Potential Changes.

Earlier, after our bioinformatic analysis of single-nucleotide polymorphisms of TATA-binding protein-binding sites within gene promoters on the human Y chromosome, we suggested that human reproductive potential diminishes during self-domestication. Here, we implemented bioinformatics models of human diseases using animal in vivo genome-wide RNA-Seq data to compare the effect of co-directed changes in the expression of orthologous genes on human reproductive potential and during the divergence of domestic and wild animals from their nearest common ancestor (NCA). For example, serotonin receptor 3A (HTR3A) deficiency contributes to sudden death in pregnancy, consistently with Htr3a underexpression in guinea pigs (Cavia porcellus) during their divergence from their NCA with cavy (C. aperea). Overall, 25 and three differentially expressed genes (hereinafter, DEGs) in domestic animals versus 11 and 17 DEGs in wild animals show the direction consistent with human orthologous gene-markers of reduced and increased reproductive potential. This indicates a reliable association between DEGs in domestic animals and human orthologous genes reducing reproductive potential (Pearson's χ2 test p < 0.001, Fisher's exact test p < 0.05, binomial distribution p < 0.0001), whereas DEGs in wild animals uniformly match human orthologous genes decreasing and increasing human reproductive potential (p > 0.1; binomial distribution), thus enforcing the norm (wild type).

Genes associated with cognitive performance in the Morris water maze: an RNA-seq study.

Learning and memory are among higher-order cognitive functions that are based on numerous molecular processes including changes in the expression of genes. To identify genes associated with learning and memory formation, here, we used the RNA-seq (high-throughput mRNA sequencing) technology to compare hippocampal transcriptomes between mice with high and low Morris water maze (MWM) cognitive performance. We identified 88 differentially expressed genes (DEGs) and 24 differentially alternatively spliced transcripts between the high- and low-MWM-performance mice. Although the sets of DEGs and differentially alternatively spliced transcripts did not overlap, both were found to be enriched with genes related to the same type of biological processes: trans-synaptic signaling, cognition, and glutamatergic transmission. These findings were supported by the results of weighted-gene co-expression network analysis (WGCNA) revealing the enrichment of MWM-cognitive-performance-correlating gene modules with very similar Gene Ontology terms. High-MWM-performance mice manifested mostly higher expression of the genes associated with glutamatergic transmission and long-term potentiation implementation, which are processes necessary for memory acquisition and consolidation. In this set, there were genes participating in the regulation of trans-synaptic signaling, primarily AMPA receptor signaling (Nrn1, Nptx1, Homer3, Prkce, Napa, Camk2b, Syt7, and Nrgn) and calcium turnover (Hpca, Caln1, Orai2, Cpne4, and Cpne9). In high-MWM-performance mice, we also demonstrated significant upregulation of the "flip" splice variant of Gria1 and Gria2 transcripts encoding subunits of AMPA receptor. Altogether, our data helped to identify specific genes in the hippocampus that are associated with learning and long-term memory. We hypothesized that the differences in MWM cognitive performance between the mouse groups are linked with increased long-term potentiation, which is mainly mediated by increased glutamatergic transmission, primarily AMPA receptor signaling.

Disruptive natural selection by male reproductive potential prevents underexpression of protein-coding genes on the human Y chromosome as a self-domestication syndrome.

BACKGROUND: In population ecology, the concept of reproductive potential denotes the most vital indicator of chances to produce and sustain a healthy descendant until his/her reproductive maturity under the best conditions. This concept links quality of life and longevity of an individual with disease susceptibilities encoded by his/her genome. Female reproductive potential has been investigated deeply, widely, and comprehensively in the past, but the male one has not received an equal amount of attention. Therefore, here we focused on the human Y chromosome and found candidate single-nucleotide polymorphism (SNP) markers of male reproductive potential. RESULTS: Examining in silico (i.e., using our earlier created Web-service SNP_TATA_Z-tester) all 1206 unannotated SNPs within 70 bp proximal promoters of all 63 Y-linked genes, we found 261 possible male-reproductive-potential SNP markers that can significantly alter the binding affinity of TATA-binding protein (TBP) for these promoters. Among them, there are candidate SNP markers of spermatogenesis disorders (e.g., rs1402972626), pediatric cancer (e.g., rs1483581212) as well as male anxiety damaging family relationships and mother's and children's health (e.g., rs187456378). First of all, we selectively verified in vitro both absolute and relative values of the analyzed TBP-promoter affinity, whose Pearson's coefficients of correlation between predicted and measured values were r = 0.84 (significance p <  0.025) and r = 0.98 (p <  0.025), respectively. Next, we found that there are twofold fewer candidate SNP markers decreasing TBP-promoter affinity relative to those increasing it, whereas in the genome-wide norm, SNP-induced damage to TBP-promoter complexes is fourfold more frequent than SNP-induced improvement (p <  0.05, binomial distribution). This means natural selection against underexpression of these genes. Meanwhile, the numbers of candidate SNP markers of an increase and decrease in male reproductive potential were indistinguishably equal to each other (p <  0.05) as if male self-domestication could have happened, with its experimentally known disruptive natural selection. Because there is still not enough scientific evidence that this could have happened, we discuss the human diseases associated with candidate SNP markers of male reproductive potential that may correspond to domestication-related disorders in pets. CONCLUSIONS: Overall, our findings seem to support a self-domestication syndrome with disruptive natural selection by male reproductive potential preventing Y-linked underexpression of a protein.

Asymmetric Conservation within Pairs of Co-Occurred Motifs Mediates Weak Direct Binding of Transcription Factors in ChIP-Seq Data.

(1) Background: Transcription factors (TFs) are main regulators of eukaryotic gene expression. The cooperative binding to genomic DNA of at least two TFs is the widespread mechanism of transcription regulation. Cooperating TFs can be revealed through the analysis of co-occurrence of their motifs. (2) Methods: We applied the motifs co-occurrence tool (MCOT) that predicted pairs of spaced or overlapped motifs (composite elements, CEs) for a single ChIP-seq dataset. We improved MCOT capability for the prediction of asymmetric CEs with one of the participating motifs possessing higher conservation than another does. (3) Results: Analysis of 119 ChIP-seq datasets for 45 human TFs revealed that almost for all families of TFs the co-occurrence with an overlap between motifs of target TFs and more conserved partner motifs was significantly higher than that for less conserved partner motifs. The asymmetry toward partner TFs was the most clear for partner motifs of TFs from the ETS (E26 Transformation Specific) family. (4) Conclusion: Co-occurrence with an overlap of less conserved motif of a target TF and more conserved motifs of partner TFs explained a substantial portion of ChIP-seq data lacking conserved motifs of target TFs. Among other TF families, conservative motifs of TFs from ETS family were the most prone to mediate interaction of target TFs with its weak motifs in ChIP-seq.

Candidate SNP Markers of Atherogenesis Significantly Shifting the Affinity of TATA-Binding Protein for Human Gene Promoters show stabilizing Natural Selection as a Sum of Neutral Drift Accelerating Atherogenesis and Directional Natural Selection Slowing It.

(1) Background: The World Health Organization (WHO) regards atherosclerosis-related myocardial infarction and stroke as the main causes of death in humans. Susceptibility to atherogenesis-associated diseases is caused by single-nucleotide polymorphisms (SNPs). (2) Methods: Using our previously developed public web-service SNP_TATA_Comparator, we estimated statistical significance of the SNP-caused alterations in TATA-binding protein (TBP) binding affinity for 70 bp proximal promoter regions of the human genes clinically associated with diseases syntonic or dystonic with atherogenesis. Additionally, we did the same for several genes related to the maintenance of mitochondrial genome integrity, according to present-day active research aimed at retarding atherogenesis. (3) Results: In dbSNP, we found 1186 SNPs altering such affinity to the same extent as clinical SNP markers do (as estimated). Particularly, clinical SNP marker rs2276109 can prevent autoimmune diseases via reduced TBP affinity for the human MMP12 gene promoter and therefore macrophage elastase deficiency, which is a well-known physiological marker of accelerated atherogenesis that could be retarded nutritionally using dairy fermented by lactobacilli. (4) Conclusions: Our results uncovered SNPs near clinical SNP markers as the basis of neutral drift accelerating atherogenesis and SNPs of genes encoding proteins related to mitochondrial genome integrity and microRNA genes associated with instability of the atherosclerotic plaque as a basis of directional natural selection slowing atherogenesis. Their sum may be stabilizing the natural selection that sets the normal level of atherogenesis.

A single ChIP-seq dataset is sufficient for comprehensive analysis of motifs co-occurrence with MCOT package.

Recognition of composite elements consisting of two transcription factor binding sites gets behind the studies of tissue-, stage- and condition-specific transcription. Genome-wide data on transcription factor binding generated with ChIP-seq method facilitate an identification of composite elements, but the existing bioinformatics tools either require ChIP-seq datasets for both partner transcription factors, or omit composite elements with motifs overlapping. Here we present an universal Motifs Co-Occurrence Tool (MCOT) that retrieves maximum information about overrepresented composite elements from a single ChIP-seq dataset. This includes homo- and heterotypic composite elements of four mutual orientations of motifs, separated with a spacer or overlapping, even if recognition of motifs within composite element requires various stringencies. Analysis of 52 ChIP-seq datasets for 18 human transcription factors confirmed that for over 60% of analyzed datasets and transcription factors predicted co-occurrence of motifs implied experimentally proven protein-protein interaction of respecting transcription factors. Analysis of 164 ChIP-seq datasets for 57 mammalian transcription factors showed that abundance of predicted composite elements with an overlap of motifs compared to those with a spacer more than doubled; and they had 1.5-fold increase of asymmetrical pairs of motifs with one more conservative 'leading' motif and another one 'guided'.

A Rat Model of Human Behavior Provides Evidence of Natural Selection Against Underexpression of Aggressiveness-Related Genes in Humans.

Aggressiveness is a hereditary behavioral pattern that forms a social hierarchy and affects the individual social rank and accordingly quality and duration of life. Thus, genome-wide studies of human aggressiveness are important. Nonetheless, the aggressiveness-related genome-wide studies have been conducted on animals rather than humans. Recently, in our genome-wide study, we uncovered natural selection against underexpression of human aggressiveness-related genes and proved it using F1 hybrid mice. Simultaneously, this natural selection equally supports two opposing traits in humans (dominance and subordination) as if self-domestication could have happened with its disruptive natural selection. Because there is still not enough scientific evidence that this could happen, here, we verified this natural selection pattern using quantitative PCR and two outbred rat lines (70 generations of artificial selection for aggressiveness or tameness, hereinafter: domestication). We chose seven genes-Cacna2d3, Gad2, Gria2, Mapk1, Nos1, Pomc, and Syn1-over- or underexpression of which corresponds to aggressive or domesticated behavior (in humans or mice) that has the same direction as natural selection. Comparing aggressive male rats with domesticated ones, we found that these genes are overexpressed statistically significantly in the hypothalamus (as a universal behavior regulator), not in the periaqueductal gray, where there was no aggressiveness-related expression of the genes in males. Database STRING showed statistically significant associations of the human genes homologous to these rat genes with long-term depression, circadian entrainment, Alzheimer's disease, and the central nervous system disorders during chronic IL-6 overexpression. This finding more likely supports positive perspectives of further studies on self-domestication syndromes.

Rare Earth Elements Alter Redox Balance in Methylomicrobium alcaliphilum 20ZR.

Background: Rare Earth Elements (REEs) control methanol utilization in both methane- and methanol-utilizing microbes. It has been established that the addition of REEs leads to the transcriptional repression of MxaFI-MeDH [a two-subunit methanol dehydrogenase (MeDH), calcium-dependent] and the activation of XoxF-MeDH (a one-subunit MeDH, lanthanum-dependent). Both enzymes are pyrroquinoline quinone-dependent alcohol dehydrogenases and show significant homology; however, they display different kinetic properties and substrate specificities. This study investigates the impact of the MxaFI to XoxF switch on the behavior of metabolic networks at a global scale. Results: In this study we investigated the steady-state growth of Methylomicrobium alcaliphilum 20ZR in media containing calcium (Ca) or lanthanum (La, a REE element). We found that cells supplemented with La show a higher growth rate compared to Ca-cultures; however, the efficiency of carbon conversion, estimated as biomass yield, is higher in cells grown with Ca. Three complementary global-omics approaches-RNA-seq transcriptomics, proteomics, and metabolomics-were applied to investigate the mechanisms of improved growth vs. carbon conversion. Cells grown with La showed the transcriptional activation of the xoxF gene, a homolog of the formaldehyde-activating enzyme (fae2), a putative transporter, genes for hemin-transport proteins, and nitrate reductase. In contrast, genes for mxaFI and associated cytochrome (mxaG) expression were downregulated. Proteomic profiling suggested additional adjustments of the metabolic network at the protein level, including carbon assimilation pathways, electron transport systems, and the tricarboxylic acid (TCA) cycle. Discord between gene expression and protein abundance changes points toward the possibility of post-transcriptional control of the related systems including key enzymes of the TCA cycle and a set of electron-transport carriers. Metabolomic data followed proteomics and showed the reduction of the ribulose-monophosphate (RuMP) pathway intermediates and the increase of the TCA cycle metabolites. Conclusion: Cells exposed to REEs display higher rates of growth but have lower carbon conversion efficiency compared to cells supplemented with Ca. The most plausible explanation for these physiological changes is an increased conversion of methanol into formate by XoxF-MeDH, which further stimulates methane oxidation but limits both the supply of reducing power and flux of formaldehyde into the RuMP pathway.

The Interplay of Chromatin Landscape and DNA-Binding Context Suggests Distinct Modes of EIN3 Regulation in Arabidopsis thaliana.

The plant hormone ethylene regulates numerous developmental processes and stress responses. Ethylene signaling proceeds via a linear pathway, which activates transcription factor (TF) EIN3, a primary transcriptional regulator of ethylene response. EIN3 influences gene expression upon binding to a specific sequence in gene promoters. This interaction, however, might be considerably affected by additional co-factors. In this work, we perform whole genome bioinformatics study to identify the impact of epigenetic factors in EIN3 functioning. The analysis of publicly available ChIP-Seq data on EIN3 binding in Arabidopsis thaliana showed bimodality of distribution of EIN3 binding regions (EBRs) in gene promoters. Besides a sharp peak in close proximity to transcription start site, which is a common binding region for a wide variety of TFs, we found an additional extended peak in the distal promoter region. We characterized all EBRs with respect to the epigenetic status appealing to previously published genome-wide map of nine chromatin states in A. thaliana. We found that the implicit distal peak was associated with a specific chromatin state (referred to as chromatin state 4 in the primary source), which was just poorly represented in the pronounced proximal peak. Intriguingly, EBRs corresponding to this chromatin state 4 were significantly associated with ethylene response, unlike the others representing the overwhelming majority of EBRs related to the explicit proximal peak. Moreover, we found that specific EIN3 binding sequences predicted with previously described model were enriched in the EBRs mapped to the chromatin state 4, but not to the rest ones. These results allow us to conclude that the interplay of genetic and epigenetic factors might cause the distinct modes of EIN3 regulation.